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A  STUDY  ON  THE  PREDICTION  OF  THE  STARTING  TRANSIENT 
PROCESS  OF  SOLID  PROPELLANT  ROCKET  MOTORS 


Chien  Tse-ch'un,  Ho  Hung-Ch'ing 


ABSTRACT 


The  interior  ballistic  P(x,  t)  model  of  the  starting  transient  process  of  the  solid  pro¬ 
pellant  rocket  motor  has  been  unproved  in  the  present  paper,  with  the  aim  of  ensuring 
better  accordance  with  the  actual  working  conditions  in  the  motor.  The  numerical  solu¬ 
tion  of  the  governing  equations,  pertaining  to  one-dimensional  unsteady  gas  dynamics, 
ntiliTfc  an  implicit  finite-difference  scheme  combined  with  the  method  of  characteristics. 
The  computed  pressure-time  curve  is  in  fair  agreement  with  the  experimentally  measured 
one.  It  is  evident  from  the  results  of  computation  that  erosive  burning  is^  main  factor 
affecting  the  starting  pressure-peak  for  a  given  propellant  geometry,  and  the  heat  exchange 
between  the  combustion-gas  and  the  solid  propellant  casts  a  great  influence  upon  the 
processes  of  ignition-induction  and  flame-spreading.  Moreover  the  rate  of  pressure  change 
also  plays  an  important  role  in  the  phases  of  flame-spreading  and  chamber-filling,  and 
should  hence  be  accounted  for  in  the  prediction.  The  computation  program  presented 
should  provide  a  general  means  for  accurate  prediction  of  the  interior  ballistic  curve  of 
starting  transient  process,  as  well  as  for  the  investigation  of  the  influences  on  it  of  the 
various  factors  involved.  t 


I .  Introduction 


The  interior  ballistic  curve  of  a  solid  propellant  rocket 
consists  of  the  starting  transient  stage,  the  operating  stage  and 
the  trailing  stage  (Figure  1) .  The  starting  transient  stage,  in 
turn,  is  made  up  of  the  ignition-induction  stage,  the  flame¬ 
spreading  stage  and  the  chamber- filling  stage.  Ignition  delay, 
rate  of  pressure  increase  and  peak  value  of  the  initial  pressure 
may  be  determined  through  computation  of  the  starting  process . 

The  purpose  of  in-depth  studies  of  the  starting  transient  process 
is:  1.  to  obtain  an  accurate  prediction  method  of  the  starting 
transient  process  to  serve  as  a  basis  for  the  choice  by  the 
designer  of  appropriate  propellant  geometry,  initial  combustion 
surface  and  shell  body  thickness,  thus  raising  as  much  as  possible 
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the  mass  ratio  of  the  motor,  reducing  the  number  of  necessary 
ground  tests,  and  thus  slashing  cost  and  shortening  the  period 
of  fabrication.  2.  to  study  the  effect  of  various  parameters 
so  as  to  provide  a  basis  for  the  control  of  the  starting  transient 
process . 


t 


Figure  1.  Schematic  diagram  of  the  operation  of  a  solid 
propellant  rocket  motor. 

1 — starting  transient  stage;  2 — operating  stage;  3 — trailing  stage; 
4 — erosive  burning;  5 — flame-spreading  stage;  6 — ignition-induction 
stage;  7 — chamber -filling  stage;  8 — peak  of  initial  pressure 

Two  models  were  adopted  abroad  in  the  60's  for  the  study  of 
the  starting  transient  process.  The  first  is  a  constant,  steady, 
one-dimensional  flow  model  (the  p(x)  model)  which  assumes  that  the 
combustion  gas  parameters,  p,  T  and  u,  vary  only  with  the  axial 
position  at  a  given  instant  of  time.  The  second  model  assumes  that 
p,  T  and  u  do  not  vary  along  the  passage,  but  only  vary  with  time 
(the  p(t)  model).  In  the  70’s,  Kuo,  K.  K.,  first  took  into  consid¬ 
eration  both  the  space  and  time  dependence  of  these  parameters 
(the  p(x,t)  model)  in  the  solution  of  the  one-dimensional  unsteady 
gas  dynamics  equations,  and  accomplished  the  computation  of  the 
starting  transient  process  of  granular  propellants  [1] .  Later, 


this  model  was  generalized  by  Pertz,  et  al.,  [2]  and  applied  to 
solid  propellant  rocket  motors.  As  the  p(x,t)  model  can  be  used 
to  compute  simultaneously  the  variation  with  time  and  space  of 
the  combustion  gas  parameters,  it  gives  a  more  nearly  true  reflec¬ 
tion  of  the  transient  processes  of  ignition-induction,  flame¬ 
spreading  and  combustion  chamber-filling.  Moreover,  in  contrast 
with  the  p(t)  model,  the  p(x,t)  model  allows  the  cross-sectional 
area  of  the  solid  propellant  passage  to  vary  along  the  axial 
direction . 

Up  until  today,  on  the  p(x)  model  and  the  p(t)  model  have  been 
used  in  this  country.  To  better  understand  the  starting  transient  1,2 
process,  we  have  improved  the  p(x,t)  model,  written  a  FORTRAN  pro¬ 
gram,  and  given  the  computed  results  for  the  cases  where  certain 
parameters  have  been  changed. 

Notations 

2 

A.  combustion  surface,  cm 

°  2 
Ap  cross-sectional  area  of  passage,  cm 

At  area  at  the  throat  of  the  engine  nozzle,  cm2 

a  coefficient  of  combustion  rate  in  the  formula  for 

combustion  rate  r  =  apn 

b  combustion  perimeter,  cm 

hc  local  convection  heat  exchange  coefficient,  kcal/cm  *sec*k 

Iu,Ip,IT  inhomogeneous  terms  of  the  governing  equations,  Eq.  (1) 

Jc  mechanical  equivalent  of  heat,  0.4267xl05  kg*cm/kcal 

ifLg  igniter  gas  mass  flow  rate,  kg* sec/cm 
n  pressure  index 

2 

p  static  pressure,  kg/cm 

2 

p.  pressure  at  the  bow,  kg/cm 

2 

q^  heat  flow  rate,  kcal/cm  *sec 

R  gas  constant  of  the  combustion  gas,  kg • cm/kg »k 

r  combustion  rate  of  solid  propellant,  cm/sec 

rQ  steady-state  combustion  rate,  cm/sec 

T  static  temperature  of  the  combustion  gas,  K 


temperature  of  the  adiabatic  flame,  K 
Tp^  initial  temperature  of  the  propellant,  K 
Tps  local  surface  temperature  of  the  propellant,  K 
Tgig  critical  ignition  temperature  of  propellant,  K 
t  time,  sec 

u  velocity  of  gas,  cm/sec 

u  .  limiting  velocity  of  erosive  burning,  m/sec 

tv  , 

Ves  volume  of  passage  upstream  xq,  cm 

x  the  axial  coordinate  of  the  point  first  ignited,  cm 

x  the  axial  coordinate  of  the  front  section  of  the  solid 

P 

propellant,  cm 

2 

a  coefficient  of  heat  diffusion,  cm  /sec 

X  heat  conduction  coefficient  (in  the  absence  of  a  sub¬ 

script,  it  refers  to  the  gas),  kcal/cm»sec*k 

p  density  (in  the  absence  of  a  subscript,  it  refers  to  the 

2  4 

gas) ,  kg  sec  cm 
Subscripts: 

e  nozzle  entrance 

es  upstream  xQ 

fr  flame  front 

i  initial  value 

ig  igniter 

II .  Physical  model 

The  physical  process  of  the  starting  transient  process  can 
be  described  as  follows: 

1.  Solid  propellant  ignition-induction  stage 

The  high  temperature  combustion  gas  from  the  igniter  pro¬ 
duces  changes  in  the  parameters  of  the  flow  field  in  the  passage, 
and  undergoes  heat  exchange  with  the  surface  of  the  solid  propellant 
When  the  net  energy  acquired  at  a  certain  location  on  the  propellant 


surface  becomes  high  enough  to  maintain  spontaneous  combustion. 


ignition  occurs  at  that  location.  It  is  usually  thought  that  the 
heat  flow  is  greatest  where  the  flame  from  the  igniter  comes  in 
contact  with  the  surface  of  the  propellant  (Figure  2,  xq) ,  and  this 
is  where  the  propellant  becomes  ignited  first.  The  position  of  xq 
may  be  computed  with  the  help  of  high  speed  photography  of  the 
shape  of  the  flame  from  the  igniter  and  knowledge  about  the  pro¬ 
pellant  geometry.  It  may  also  be  obtained  from  measurements  of  the 
heat  flow  in  the  ignition  process. 


raaasasssa  Manas 


fmmx 


Figure  2.  Schematic  diagram  of  a  solid  propellant 
rocket  motor 


2 .  Flame  spreading  stage 

After  the  propellant  is  ignited  locally  on  the  surface,  the 
additional  increase  in  mass  and  energy  from  the  propellant  enhances 
the  transfer  of  energy  toward  the  unignited  parts  of  the  surface, 
causing  the  surface  of  the  propellant  close  to  the  flame  tip  to  be 
ignited  rapidly.  Thus,  the  flame  spreads  from  xq  upstream  and  down 
stream,  until  the  entire  surface  of  the  propellant  is  ignited. 

3.  Chamber-filling  stage 


After  the  entire  surface  of  the  propellant  is  ignited,  erosive 
burning  and  the  instantaneous  rate  of  change  in  pressure  cause  the 
rate  of  formation  of  the  combustion  gas  to  increase  continually. 
Pressure  in  the  combustion  chamber  rapidly  builds  up  and  an  initial 
pressure  peak  is  produced. 


In  summary,  the  physical  factors  we  have  taken  into  consider¬ 
ation  include:  (1)  positioning  of  the  igniter  in  the  propellant 
passage  at  the  bow  of  the  motor,  (2)  the  heat  exchange  process  and 
friction  between  the  combustion  gas  and  the  combustion  surface  of 
the  solid  propellant,  (3)  ignition  of  the  propellant  with  the 
flame  spreading  upstream  and  downstream,  (4)  the  effect  of  the  ero¬ 
sive  burning,  (5)  the  effect  of  the  instantaneous  rate  of  change 
of  pressure  and  (6)  the  cross-sectional  area  of  the  passage  being 
continually  variable  along  the  axial  direction.  Among  these.  Cl)  , 
(3)  and  (6)  constitute  the  improvements  we  have  made  in  this  paper. 

In  this  manner,  the  physical  model  considered  is  a  fairly 
close  representation  of  the  actual  operating  conditions  of  the 
motor. 


Ill .  Analytical  model 
1 .  The  governing  equations 

Take  the  point  xq,  where  the  flame  from  the  igniter  impinges 
on  the  surface  of  the  propellant,  and  the  nozzle  entrance  xg  to  be, 
respectively,  the  left  and  right  boundary  of  the  flow  field  (Figure 
2).  Assume  that:  Cl)  the  effect  of  the  solid  phase  on  the  flow  is 
neglected;  (2)  the  combustion  gas  obeys  the  ideal  gas  law,  and  the 
isobar ic  specific  heat,  specific  heat  ratio,  average  molecular 
weight  and  gas  constant  are  the  same  for  the  combustion  gases  of 
the  igniter  and  the  propellant,  and  remain  constant;  (3)  combustion 
of  the  propellant  is  limited  to  the  thin  layer  of  gaseous  phase 
close  to  the  surface  of  the  propellant,  and  the  addition  of  combus¬ 
tion  gas  to  the  passage  has  zero  axial  velocity  component.  After 
listing  and  transforming  the  equation  of  continuity,  the  momentum 
conservation  equation  and  the  energy  conservation  equation,  we 
obtain  the  following  set  of  governing  equations  for  one-dimensional 
unsteady  gas  dynamics. 


5+*?- £i*s£L  .. 

&  dx^B  p  9x  A,  Jur~A^fc~2mI‘ 

%  +  rp  |r+Mff - ^  dix-+(r-\)(jcgCpT,+^) 

‘  ^rr-(r_1)^7(/c«'_2^ir 

•  ['**P#r  (  r,  +  27^J  - f  )r - (r- 1 )(/ Al - ^r^/cH.) J./r 


The  initial  conditions  are 

p-pt*  *-«(,  r-r, 


(2) 


Equation  (1)  is  a  set  of  pseudo-linear,  inhomogeneous,  hyper¬ 
bolic  partial  differential  equations.  Its  three  characteristic 
directions  and  the  corresponding  compatible  equations  are,  respect 
ively,  the  characteristic  line  I  along  the  right  side 


[*i— * 

the  characteristic  line  II  along  the  left  side 

[fl— c 

the  locus  line  III  along  the  fluid 

Tdx- 1  r dTl  V-l  T  dp  ,  T  V—  1  T  . 

L*J.u“  V  p  dt+Ir~  V  pl» 


The  equation  for  the  rate  of  combustion  is  taken  to  be  the 
erosion  relation  [3]  proposed  by  Kuo,  K.  K.,  on  the  basis  of  a 
thermo-chemical  model  of  gas  dynamics.  We  also  take  into  consider 
ation  the  rapid  rise  of  pressure  in  the  combustion  chamber  during 
chamber-filling.  Combining  the  above  with  the  effect  of  rate  of 
pressure  variation  dp/dt  on  the  rate  of  combustion  [4],  we  obtain 

r-r,[l+JClNi-«,ryr.+^J]  (4) 

The  temperature  of  the  surface  of  the  propellant  is  computed 
from  the  following  ordinary  differential  equation: 


dT,Jdt-ia,M:a'-Tt'y/m',r<Jt.-TpiX2T-T,.-Ttln  (5) 

The  initial  condition  is  T  (x,0)  =  T  ..  When  T  £  T  .  on  the 

ps  pi  ps  sig 

computed  cross-section,  the  propellant  is  ignited  at  that  position. 


Boundary  conditions 


Take  the  section  x^  %  x  $  xQ  to  the  left  of  the  left  boundary 
xq  as  the  front  section  (Figure  2) .  The  flow  field  is  rather  com¬ 
plex  there.  Assume  that  the  front  section  has  p  uniform  distribu¬ 
tion  of  pressure  and  temperature.  Computation.*  iave  shown  that 
variations  in  temperature  and  pressure  in  the  :  „nt  section  are  not 
large.  We  obtain  from  the  continuity  equation  id  the  energy  con¬ 
servation  equation  applied  to  the  front  sectioi 


dpJdt-[r*l,(t')YgRTi,.+rA>..p„YgRTf 

A  (V— 1  )Apx»  s  ~\/v 
-Yp..»*.Ap*» - 2gRT,~  Pu  _]/  " 

dTJdi  =[  gRTJp'.thi,  (VT 

•  gRT^VT,- Tm, ) -  (  r  - 1  %W-.(  T ». +|^)]/K" 


(6) 


The  initial  conditions  are 

These  plus  the  equations  corresponding  to  the  characteristic  line 
along  the  left  side  in  equation  (3)  constitute  the  left  boundary 
conditions.  denotes  that  part  of  the  surface  (6)  in  the  front  sec¬ 

tion  that  is  already  undergoing  combustion.  It  is  found  from  the 
distance  by  which  the  flame  has  spread  to  the  left. 

The  right  boundary  conditions  are  formed  by  taking  the  equa¬ 
tions  for  the  one-dimensional,  isentropic,  constant  steady  choked 
or  unchoked  (before  the  nozzle  becomes  choked)  flow  of  the  nozzle, 
combined  with  the  equations  in  (3)  that  correspond  to  the 
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characteristic  line  along  the  right  side  and  those  corresponding 

to  the  locus  line  of  the  fluid.  Before  the  choke  lid  is  opened, 

=  0. 

e 

IV.  Numerical  solution 


An  implicit  finite-difference  scheme  is  combined  with  the 
method  of  characteristics  to  solve  the  partial  differential  equa¬ 
tions  whose  inhomogeneous  terms  have  been  linearized.  We  finally 
arrive  at  a  set  of  linear  algebraic  equations  which  enable  one  to 
solve  for  the  gas  flow  parameters  at  the  various  nodes  at  time 
j  +  1  given  the  parameters  for  the  various  nodes  at  the  time  j . 


(A,  C, 

Bt  At  C, 

B.  A.  C. 

iZx 


Fx-BV!" 

Ft 

K 

fZ-cZ 


(7) 


In  the  above,  U^+1  and  U^+1  are,  respectively,  the  combustion  gas 
parameters,  u,  p  and  T  on  the  left  and  right  boundary  at  time  j  +  1, 
and  are  solved  from  the  left  and  right  boundary  conditions,  respect¬ 


ively.  A,  B  and  C  are  3x3  submatrices,  and  U  and  F  are  3x1  sub¬ 


matrices  . 


The  left  and  right  boundary  conditions  have  also  been  trans¬ 
formed  into  algebraic  equations  in  a  similar  manner. 

Two  points  need  to  be  mentioned  here  about  the  procedure  of 
solving  these  equations: 

1.  After  the  propellant  first  becomes  ignited  at  xq,  the 
flame  spreads  simultaneously  upstream  and  downstream.  One  may  pro¬ 
ceed  with  the  computation  as  follows.  First,  assume  an  acceleration 
for  the  spreading  of  the  flame  toward  the  left.  Use  this  to  cal¬ 
culate  the  area  A^  of  the  front  section  that  has  already  been 
ignited.  Substitute  this  into  the  left  boundary  condition  which 
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then  is  used  along  with  the  right  boundary  condition  to  solve  the 
computational  equations,  Equation  (7),  for  the  governing  equations. 

After  the  acceleration  of  the  spreading  of  the  flame  towards  the 
right  from  point  x  to  a  neighboring  node  downstream  has  been  com¬ 
puted,  the  process  of  computation  is  started  anew  from  the  time  xQ  < 

is  ignited,  the  acceleration  of  flame  spreading  toward  the  right 
being  used  to  correct  that  toward  the  left.  This  procedure  is 
repeated  until  two  successive  acceleration  values  obtained  for  the 
spreading  of  flame  toward  the  left  are  very  close  to  each  other. 


2.  By  taking  into  consideration  the  increase  in  rate  of  com¬ 
bustion  due  to  the  rate  of  change  of  pressure  dp/dt  in  the  combus¬ 
tion  rate  relation  (4) ,  we  will  amplify  any  small  error  in  pressure 
in  the  computation  and  cause  sharp  oscillations  in  the  p-t  curve 
for  the  flame-spreading  stage.  This  difficulty  is  overcome  by 
adopting  large  steplength  Ap/At  instead  of  dp/dt,  and  resorting  to 
the  method  of  iterative  correction. 


/49 

V.  Results  of  computation,  analysis  and  conclusion 

The  numerical  computation  was  carried  out  for  the  G5  motor. 

See  Table  1  for  the  numerical  inputs . 


TABLE  1.  Numerical  inputs 


Geometric  para- 


Properties  of  the  propellant 


Combustion  gas 
parameters 


Xp 

-3.5 

0. 1286X10-* 

*tr 

63.90 

r 

1.165 

x% 

10.0  ; 

C,r 

0.29 

K„ 

0.5 

w 

27.30 

x r 

123.7 

i, 

1.8X10-* 

K , 

0 .2144  X 10“: 

R 

3100.0 

At 

50.72 

Of 

0.0023 

K, 

0.35 

c. 

0  $023 

V.„ 

1391.6 

T..t 

293.0 

K . 

0.69 

Pr 

0.6647 

6/ x) 

11.64  ' 

a 

0.2106 

r.„ 

$60.0 

T„. 

2250  0 

Ap ,( x) 

113.54  ; 

m 

0.1163  . 

«. 

0.01 

In  the  computations,  we  have  used  different  values  for  the 
initial  values,  the  adiabatic  flame  temperature  of  the  propellant, 
the  mass  flow  rate  for  the  igniter,  erosive  combustion  rate,  heat 
exchange  coefficient  and  ignition  temperature  of  the  propellant. 
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All  the  computations  are  steady  and  convergent.  The  results  of 
these  computations,  analysis  and  conclusions  arrived  at  are  pre¬ 
sented  below. 

1.  Figure  3  gives  the  p-t  curve  for  the  starting  transient 
process.  It  can  be  seen  that  the  curve  that  takes  into  account 
the  effect  of  dp/dt  is  closer  to  the  experimental  curve  for  the 
G5  motor.  Therefore,  the  computational  procedure  presented  in 
this  paper  can  be  used  to  accurately  predict  the  starting  transient 
process,  and  to  study  the  effects  of  the  various  factors  involved. 

In  addition,  the  curve  also  exhibits  the  phenomenon  of  sponta¬ 
neous  acceleration  of  flame-spreading. 

2.  It  can  be  seen  from  Figures  3  and  4  that  dp/dt  has  a  marked 
effect  on  the  p-t  curve  and  the  gas  flow  parameters,  and  should, 
therefore,  be  properly  taken  account  of.  It  can  be  seen  that, 
owing  to  the  effect  of  the  instantaneous  changes  in  the  pressure, 
the  velocity  of  flame-spreading  is  increased,  the  rate  of  pressure 
increase  is  raised,  the  time  required  to  reach  the  initial  pressure 
peak  is  shortened,  and  the  combustion  gas  velocity  at  every  point 
of  the  flow  field  is  increased.  This  is  because  when  the  pressure 
increases  sharply,  the  combustion  rate  of  the  propellant  is  in¬ 
creased,  the  rate  of  mass  flow  is  increased,  and  the  effect  of  ero¬ 
sive  burning  is  enhanced. 

3.  The  erosive  property  of  the  propellant  has  a  large  effect 
on  the  initial  pressure  peak  value  as  well  as  the  rate  of  pressure 
increase  (Figure  5) .  In  the  study  of  the  parameters,  it  was  found 
that  for  a  given  geometry  of  the  propellant,  erosive  burning  is 
the  major  factor  affecting  the  initial  pressure  peak  value. 

4.  Figure  6  reflects  the  strong  influence  of  the  heat  exchange 
between  the  combustion  gas  and  the  unignited  propellant  surface  on 
flame-spreading  and  the  rate  of  pressure  increase.  The  higher  the 
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Figure  3.  Comparison  of  computed  p-t  curves  with  experimental 
pt-curve 

1 — computed  curves;  2 — experimental  curve;  3 — effect  of  dp/dt  taken 
into  account;  Q  — peak  value;  4 — effect  of  dp/dt  not  taken  into 
account;  A — peak  value;  5 — first  point  ignited;  6--second  point 
ignited;  7 — third  point  ignited;  8 — fourth  point  ignited;  9 — end 
of  the  flame-spreading  stage;  10 — eleventh  point  ignited  (m=l0) 
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Figure  4.  Distribution  of  velocity  along  the  passage  at  various 
moments 

1 — effect  of  dp/dt  is  not  taken  into  account.  t=0.0705;  2--flame 
front  has  reached  this  point;  3 — effect  of  dp/dt  taken  into  account 
4 — effect  of  dp/dt  not  taken  into  account 
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Figure  5.  Effect  of  erosive  Figure  6.  Effect  of  heat  exchange 
property  of  propellant  on  coefficient  on  flame-spreading 

the  starting  transient 

process  1 — end  of  flame-spreading  stage 


rate  at  which  energy  is  acquired  by  the  propellant,  the  faster 
the  flame -spreading,  and  the  greater  the  rate  of  pressure  increase 
Hence,  simulated  tests  of  the  starting  transient  process  should 
also  include  a  simulation  of  the  mode  of  flame-ejection  of  the 
igniter,  rate  of  mass  flow,  combustion  properties  of  the  fuel  in 
the  igniter,  etc. 
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